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ABSTRACT 

We present a Chandra ACIS observations of N157B, a young supernova remnant (SNR) located in 
the 30 Doradus star-formation region of the Large Magellanic Cloud. This remnant contains the most 
energetic pulsar known (J053747.39-691020.2; = 4.8 x 10^* ergs s"^), which is surrounded by a 
X-ray bright nonthermal nebula that likely represents a toroidal pulsar wind terminal shock observed 
edge-on. Two of the eight point-like X-ray sources detected in the observation are shown to have 
near-IR and optical counterparts (within O'.'S offsets), which are identified as massive stellar systems 
in the Cloud. We confirm the non-thermal nature of the comet-shaped X-ray emission feature and 
show that the spectral steepening of this feature away from the pulsar is quantitatively consistent 
with synchrotron cooling of shocked pulsar wind particles flowing downstream at a bulk velocity close 
to the speed of light. Around the cometary nebula we unambiguously detect a spatially- resolved 
thermal component, which accounts for about 1/3 of the total 0.5 - 10 keV flux from the remnant. 
This thermal component is distributed among various clumps of metal-enriched plasma embedded in 
the low surface brightness X-ray-emitting diffuse gas. The relative metal enrichment pattern suggests 
that the mass of the supernova progenitor is 2QMq. A comparison of the X-ray data with Hubble 
Space Telescope optical images now suggests that the explosion site is close to a dense cloud, against 
which a reflection shock is launched. The interaction between the reflected material and the nebula 
has likely produced both its cometary shape and the surrounding thermal emission enhancement. 
SNR N157B is apparently expanding into the hot low-density interior of the surrounding superbubble 
formed by the young OB association LII99, as revealed by Spitzer mid-infrared images. This scenario 
naturally explains the exceptionally large sizes of both the thermal and nonthermal components as 
well as the lack of an outer shell of the SNR. However, the real situation in the region is likely to 
be more complicated. We find that a partially-round-shaped soft X-ray-emitting clump with distinct 
spectral properties may result from a separate oxygen-rich remnant. These results provide a rare 
glimpse into the SNR structure and evolution in a region of recent star-formation. 
Subject headings: pulsars: general — pulsars: individual (PSR J0537— 6910) — X-rays: general — 

supernova remnant: individual (N157B, SNR 0538-691) — galaxies: individual 

(Large Magellanic Cloud) 
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1. INTRODUCTION 

The evolution of a supernova remnant (SNR) depends 
sensitively on its environment. Most massive stars are 
expected to be born in OB associations and end their 
lives in an environment altered by strong ionizing radi- 
ation and mechanical energy feedback from fast stellar 
winds and supernovae (SNe). This feedback produces 
superbubbles filled with low-density hot gas. An SNR 
inside such a superbubble will be difficult to observe 
(e.g., Tang & Wang 2005). A classical example is the 
Crab SNR, which shows little detectable emission from 
its blastwave-heated gas and is believed to be expanding 
in a low-density medium. Without the synchrotron radi- 
ation from its pulsar wind nebula (PWN), the Crab SNR 
would not have been identified. Similar SNRs dominated 
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by radiation from PWNe are called "Crab- like". Phys- 
ically, Crab-like SNRs are not very different from the 
"composite" remnants, which show comparable thermal 
and nonthermal components or show a shell-like struc- 
ture with a centrally-filled morphology. The increasing 
sensitivity of X-ray observations has made it possible to 
detect more shell-like and thermal features around Crab- 
like SNRs. For example, the Crab-like SNR G21. 5-0.9 
has recently been shown to be surrounded by ejecta- or 
blastwave-dominated features (Slane et al. 2000; Mathe- 
son & Safi-Harb 2005). In this paper, we present evidence 
that the Crab-hke SNR N157B, also known as SNR 0538- 
691, exhibits thermal X-ray emission as well and is thus 
a composite SNR. 

SNR N157B is particularly interesting, because of its 
co-existence with a young OB association LII99 in the 
Large Magellanic Cloud (LMC) (Lucke & Hodge 1970; 
Chu et al. 1992). The SNR has long been classified to 
be Crab-like and contains the fastest known young pul- 
sar PSR J0537-6910 with a period of 16 ms (Marshall 
et al. 1998). The SNR appears to have an exceptionally 
large physical size in both radio and X-ray 20 — 30 pc; 
Lazendic et al. 2000; Wang & Gotthelf 1998, WG98 here- 
after), although the outer boundaries have not yet been 
well determined. Of particular interest is how the struc- 
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ture and evolution of SNR N157B is affected by the pres- 
ence of the OB association LH99. Is the SNR expanding 
in a low-density hot medium, as may be expected in a 
superbubble created by the OB association? Addressing 
such questions is important in understanding both the 
fate of the feedback from massive stars and the physics 
involved in the evolution of such an SNR (e.g., Tang & 
Wang 2005). 

Previous X-ray studies based on data with low angular 
resolutions were concentrated on the overall spectral and 
morphological characterization of N157B. Using ASCA 
and ROSAT data, WG98 found that the X-ray emission 
arises from a bright cometary-shaped nebula, headed by 
a compact source embedded in a bright bar of nonther- 
mal emission, with a long tail of diffuse emission trail- 
ing on one side aligned perpendicular to the bar. This 
morphology is reminiscent of those associated with trail- 
ing plerions in SNRs G327. 1-1.1 (Sun, Wang, & Chen 
1999), W44 (Frail et al. 1996), and IC443 (Olbert et al. 
2001; Gaensler et al. 2006), as weh as the PWN around 
PSR B1757-24 near SNR G5.4-1.2 (the Duck X-ray 
Nebula; Kaspi et al. 2001). Using a Chandra HRC obser- 
vation with high angular and timing resolutions, Wang 
et al. (2001) were able to pinpoint (to arcsec precision) 
the location of PSR J0537-6910, which is still undetected 
in any other wavelength band (e.g. Mignani et al. 2005; 
Crawford et al. 2005). The nebula was shown to con- 
sist of primarily a ^ 0.6 pc xl.7 pc compact bar-like 
feature around the pulsar and a tail 5 pc long. This 
morphology indicates that the bar represents the reverse 
shock of a toroidal wind from PSR J0537— 6910 and that 
the ambient pressure confinement of the nebula is largely 
one-sided. 

The one-sided confinement of the N157B nebula is 
thought to be a result of the pulsar motion relative to 
the ambient medium. This hypothesis seems to be con- 
sistent with the displacement of the pulsar from an ex- 
tended radio continuum emission peak, mostly represent- 
ing the accumulation of pulsar wind particles, assumed 
to be mainly electrons and positrons (WG98). Based on 
2-D hydrodynamic simulations, van der Swaluw, Downes, 
& Keegan (2004) and van der Swaluw (2004) have con- 
cluded that a typical PWN evolves through three major 
stages: 1) expanding supersonically in freely expanding 
SN ejecta; 2) interacting with the reverse shock of the 
SNR, oscillating and relaxing to a subsonic expansion 
in the reverse-shock heated ejecta; and 3) heading into 
the general ISM, supersonically. In both the first and 
the third stages, one may expect a bow shock because of 
supersonic motion. They suggest that the PWN of the 
N157B nebula is currently in the second stage; the tail 
of the nebula is a pulsar wind bubble that was elongated 
owing to the pulsar motion and was later crushed by the 
reverse shock of the SNR. The reverse shock heats a shell 
of ejecta swept-up by the expansion of the PWN and 
bounces back into the already heated remnant materials, 
forming a parabolic layer (see Fig. 2 in van der Swaluw 
2004). In this scenario, the relative motion between the 
pulsar and its ambient medium does not need to be su- 
personic. Of course, the reality could be much more com- 
plicated than what is assumed in the simulations. The 
region shows clearly large density and temperature in- 
homogeneities that could have significantly affected the 
dynamics of both the SNR and the PWN (Fig. [lli,b) . In 



short, the exact nature of the confinement as well as the 
physical properties of the nebula are still very uncertain. 

Existing studies have also shown evidence for a weak 
thermal emission component of N157B. Much of the evi- 
dence comes from the spectral decomposition of the over- 
all emission from the SNR. The thermal emission arises 
over a region more extended than the PWN, but only 
very limited spatial analysis has been presented previ- 
ously (WG98; Dennerl et al. 2001; Townsley et al. 2006). 
It is thus of great interest to examine the physical, chemi- 
cal, and morphological properties of the thermal emission 
and to determine how they are related to the PWN and 
to the OB association and superbubble environment. 

To obtain a fresh look of the above issues, we have ac- 
quired an on-axis Chandra ACIS-S observation of N157B 
(Fig. 121), allowing for a detailed spatially-resolved X-ray 
spectroscopy of the SNR. We have also analyzed the X- 
ray data together with complementary archival images 
from recent Hubble Space Telescope {HST) Advanced 
Camera for Surveys (ACS) and Spitzer Space Telescope 
Infrared Array Camera (IRAC) observations. In § 2, we 
briefly describe the observation and data calibration, and 
present our analysis and results. The effects of the su- 
perbubble and density inhomogeneity, the origin of the 
thermal emission, and the nature of various nonthermal 
features of the PWN are discussed in § 3. We summa- 
rize our results and conclusions in § 4. Throughout this 
paper, we adopt the distance of N157B as 50 kpc (thus 
1" corresponds to 0.25 pc). Statistical errors are all pre- 
sented in the 90% confidence level. 

2. OBSERVATIONS AND DATA ANALYSIS 

The X-ray study reported here is based on our Chan- 
dra observation of N157B (ObsID 2783). The data were 
collected in two exposures, obtained on 2002 August 23- 
24 and on 2002 October 28-30, for a combined duration 
of 53 ks. The data from the two exposures were merged 
through the pipeline processing and no additional align- 
ment is attempted. Because of the separation in time, the 
relative roll angle of these exposures is 65deg (Fig. I^J. 
Data were collected with ACIS-S S3 in the focal plane us- 
ing the 1/4-chip sub-array mode (with a frame readout 
time of 0.8 s) to minimize photon pile-up associated with 
the bright pulsar. The first part of the observation (30 
ks exposure) has been used by Mignani et al. (2005) in a 
primarily optical study of PSR J0537— 6910, including a 
search for its optical counterpart and an X-ray spectral 
analysis. 

We reprocess the event files (from Level 1 to Level 2) 
for the observation using the updated CIAO data pro- 
cessing software. After removing fiares with count rates 
greater than 1.2 times the mean light curve value, a net 
exposure of 48 ks remains and is used for subsequent 
analysis. For the spectral analysis, individual spectra are 
adaptively binned to achieve a background-subtracted 
signal-to-noise ratio of 5, unless otherwise noted. The 
XSPEC spectral fitting software is used throughout. For 
the foreground absorption, we use the cross-sections of 
Morrison & McCammon (1983) and assume solar abun- 
dances. 

Fig-Elpresents the overall Chandra X-ray intensity im- 
age in the soft (0.3 - 1.5 keV) and hard (1.5 - 7.0 keV) 
energy bands, which clearly demonstrate the energy de- 
pendence of the SNR morphology. The images do not 
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Fig. 1. — (a) A Spitzer global view of the N157B region in mid-infrared (IR) (8 ^m; Fazio et al. 2004) with overlaid Chandra ACIS-S 
intensity contours at 2, 3, 4, 6, 9, 13, 18, 25, 50, 100, 200, 500, and 1000 xlO~^ counts s ~^ arcmin"^. The outer boundary of the contour 
map shows the region covered in (b). The small square boxes mark discrete X-ray sources detected, (b) Color-coded multi-band ACIS-S 
intensity images of N157B in the 0.3-0.7 keV (red), 0.7-1.5 keV (Green), and 1.5-7 keV (blue) bands, (c) ACIS-S intensity contours at 12, 
20, 30, 50, 100, 250, 800, 2500, 10000 xlO'^ counts s"! arcmin-^ overlaid on an HST ACS close-up of the N157B central region in optical: 
F814W (color-coded in red), F555W (green), and F435W (blue), (d) Equivalent width (EW) images of the emission lines: O (in red), Ne 
(green), and Mg (blue) in square-root brightness scales; the line energies are listed in Table [Tl and the image for each band is adaptively 
rebinned to include at least 10 counts per pixel. The X-ray images in (a) and (b) are smoothed with the CIAO routine csmootii, whereas 
the image in (c) is smoothed with a Gaussian kernel adjusted adaptively to achieve a count-to-noise ratio of 6. 



show a typical rim-brightened outer SNR sheU; there- 
fore, the overall size of the diffuse X-ray emission is un- 
certain. Here we will adopt a diameter of the remnant 
as ~ 100". The diffuse emission also exhibits various 
lumpy sub-structures. Intensity maps and correspond- 
ing exposure maps were also generated in four broad en- 
ergy bands, 0.3 - 0.7 keV, 0.7 - 1.5 keV, 1.5 - 3.0 keV, 
and 3.0 - 7.0 keV (Wang 2004). These maps are used to 
create the flat-fielded count rate images, which are typi- 
cally smoothed with an adaptive filter (using the program 
csmooth implemented in the CIAO software package) to 
achieve a broad-band (0.3 -7 keV) signal-to-noise ratio of 
- 3. 



To compare the distribution of the X-ray emission with 
its interstellar environment of N157B, we have used HST 
ACS/WFC and Spitzer IRAC images. The pipeline- 
processed ACS images from HST program 9471 (PI: 
Mignani), taken with the F435W, F555W, and F814W 
filters, were retrieved from the MAST archive and are 
presented without further processing. Spitzer observed 
the 30 Dor region with the IRAC (Fazio et al. 2004) 
as part of programs 63, 1032, and 20203 (Pis: Houck, 
Brandl, and Meixner, respectively). Each program used 
the IRAC mapping mode to obtain maps of the region 
with multiple dithered 12 s exposures at each location in 
the map. The Basic Calibrated Data (BCDs) from these 




Fig. 2. — The raw data image of the Chandra ACIS-S observation 
of SNR N157B. The coordinate axes are for R.A. {J2000) and Dec. 
(J2000), respectively. 



TABLE 1 

Energy Bounds of the EW images 



Elements 


Line 


Low^ 


High'' 




(eV) 


(eV) 


(eV) 


O 


590-740 


300-550 


1160-1260 


Ne 


870-1140 


750-870 


1160-1260 


Mg 


1290-1430 


1160-1260 


1620-1720 



^ The low- and high-energy ranges around the selected line en- 
ergies used to estimate the underlying continuum. 



programs were combined using standard routines in the 
MOPEX software package to obtain images with ^180 s 
integration time at all locations where the programs over- 
lap. More information on the instruments, pipeline pro- 
cessing, and the MOPEX package can be found at the 
Spitzer Science Center's Observer Support website^. 

Fig.fflpresents various perspectives of SNR N157B and 
its environment. The X-ray photon energy dependence 
is illustrated in tricolor images, for broad-band inten- 
sities and for equivalent widths (EWs) of key emission 
lines listed in Table 1. The comet-shaped nebula, as 
presented in WG98 and Wang et al. (2001), stands out 
as the brightest and hardest X-ray feature in the field. 
Surrounding the nebula is a large-scale diffuse emission 
with relatively low-surface brightness and soft spectral 
characteristics. These smoothed X-ray images are also 
compared with both optical and mid-IR maps to illus- 
trate their complex relationship. 



Point sources in the field were detected following the 
procedure detailed in Wang (2004) using a combination 
of algorithms: wavelet, sliding box, and maximum like- 
lihood centroid fitting in the above energy bands. The 

® http:/ /ssc. spitzer. caltech.edu/ost. 
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source count rate is estimated based on the net counts 
within the 90% energy-encircled radius (EER) of the 
telescope/ ACIS point-spread function (PSF; Jerius et al. 
2000). Table[5]lists our detected point-like sources whose 
positions are marked in Fig. Because of the exposure 
correction for the exposure is very small in the subarry 
(smaller than the systematic errors involved in calculat- 
ing the count rates), the actual count rates without cor- 
rection are included in the table. 

We search for possible counterparts in the SIM- 
BAD, NED, and Two Micron All Sky Survey (2MASS) 
databases (Skrutskie et al. 2006) as well as in various op- 
tical images. We find two optical stellar counterparts, 1- 
28 (type 05 f*, V=14.26) and 1-98 [05 V((f)), V=13.92; 
Schild & Testor 1992], for the X-ray sources # 3 and 7, 
respectively (TableEjl. The most accurate (O'.'l RMS) po- 
sitions of these two stars are obtained from 2MASS: R.A. 
(J2000), Dec. (J2000) = 5'^37"37r97, -69°10'14'.'7 for 1- 
28 and 5'>37"59?43, -69°9'l'.'l for 1-98; the offsets from 
the respective X-ray sources are O'.'l and 0'.'5, consistent 
with their position uncertainties (Table |2Jl. Therefore, 
the absolute astrometry accuracy of the X-ray observa- 
tion is likely better than 0'.'5. Except for the pulsar (# 
4), the other detected X-ray sources are all quite faint. 
Assuming the X-ray spectrum as a thermal plasma with a 
temperature of 0.6 keV (for soft case) and 2 keV (for hard 
case), respectively, typical for normal early-type stars, 
and a foreground absorption of '--^ 6 x 10^^ cm~^ 
in the region, we estimate the 0.3-8 keV luminosities 
of the two X-ray sources as ^ (1.2-3.6) x lO'^^ergss"^ 
based on the count rates in Table 13 These luminosities 
are higher than 10"^^ ergs s~^, indicating that the X-ray 
sources most probably represent massive colliding stellar 
wind binaries, while they may still possibly be individ- 
ual early 0-type or WR stars (e.g., Oskinova 2005). The 
source # 6, very soft and without an optical counterpart, 
may be a foreground cataclysmic variable. Other sources 
are at larger off-pulsar distances and are probably also 
not related to N157B. 



2.1. Spectral Properties 

We extracted spectra for analysis from the various fea- 
tures marked in the raw count image shown in Fig. ^ 
This figure shows some details that have been smoothed 
out in the previous images. The overall spectrum of the 
SNR is taken from a large rectangular region (labeled 
"ENT" in Fig. 21, including the pulsar. The backgroimd 
spectrum is extracted from the southern and northern 
rectangular regions (Bl and B2) marked in Fig. 0] Be- 
cause of the uncertainty of the actual extent of SNR, the 
background regions (Bl and B2) may still be contam- 
inated by the diffuse emission from the remnant. The 
background-subtracted spectrum is shown in Fig. |S1 

Though nearly featureless above ~ 2keV, the spec- 
trum shows evidence for the Ka emission lines of O VIII 
(- 0.65 keV), NelX (- 0.92 keV) andNeX(- 1.05 keV), 
Mg XI (~ 1.34 keV), and possibly Si XIII (~ 1.87 keV). 
We fit the spectrum with a power law plus a non- 
equilibrium ionization (NEI) thermal plasma that char- 
acterizes the parameters and ionization timescale {n^ti) 
of the hot gas (Borkowski, Lyerly, & Reynolds 2001). 
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TABLE 2 
Chandra SOURCE List 



Source 
(1) 


CXOU Name 
(2) 


(") 
(3) 


CR { cts ks-i) 
(4) 


HR 

(5) 


HRl 

(6) 


Flag 
(7) 


1 


J053715. 18-691316.6 


0.7 


2.49 ± 0.40 


0.67 ±0.13 




B 


2 


J053727.56-691218.4 


0.6 


1.22 ±0.27 






B 


3 


J053737.97-691014.8 


0.4 


0.32 ±0.10 






B 


4 


J053747.39-691020.2 


0.2 


205.83 ±2.30 


0.36 ±0.01 


0.92 ±0.01 


H 


5 


J053751. 79-691143.3 


0.7 


0.30 ±0.11 






H 


6 


J053759. 15-691018.4 


0.3 


0.83 ±0.15 


-0.89 ±0.13 


0.97 ±0.06 


B 


7 


J053759. 48-690901. 6 


0.4 


0.62 ±0.20 






B 


8 


J053801. 42-691229.1 


0.6 


0.74 ±0.18 






B 



Note. — The definition of the bands: 0.3-0.7 (SI), 0.7-1.5 (S2), 1.5-3 (HI), and 3-7 keV (H2). In addition, S=S1+S2, H=H1+H2, and 
B=S+H. Column (1): Generic source number. (2): Chandra X-ray Observatory (unregistered) source name, following the Chandra naming 
convention and the lAU Recommendation for Nomenclature (e.g., http://cdsweb.u-strasbg.fr/iau-spec.htmll. (3): Position uncertainty, 
including an la statistical error calculated from the maximum likeUhood centroiding and an approximate off-axis angle (r) dependent 
systematic error 0'.'2 + l''4(r/8')'^ (an approximation to Fig. 4 of Feigelson et al. 2002), which are added in quadrature. (4): On-axis 
source broad-band count rate — the sum of the exposure-corrected count rates in the four bands. (5-6): The hardness ratios defined as 
HR = (H - S2)/(H + S2), and HRl = (82 - Sl)/S, hsted only for values with uncertainties less than 0.2. (7): The labels "B", "S", or "H" 
mark the bands in which a source is detected with the most accurate position that is adopted in Column (2). 



Our fit allows the abundances of the above significant 
line-emitting elements to vary, while other metal ele- 
ments are set to be 0.3 solar, typical for the LMC (Russel 
& Dopita 1992). 

The fit results are summarized in Table The lumi- 
nosity of the nonthermal 0.5-10 keV emission, including 
the contributions from PSR J0537-6910 and the PWN, 
is Lnt ~ 4.4 X 10'^^ ergs s~^, twice that of the thermal 
emission Lth ~ 2.2 x 10^® ergs s~^. Clearly, the overall 
spectrum of N157B is dominated by the power-law emis- 
sion, particularly at high energies (> 2 keV). As a com- 
parison, the overall luminosity of N157B is independently 
estimated as 4.2 x 10"^^ ergs s^^ (0.5 - 8 keV) from the 
Chandra ACIS-I observation applying a coUisional ion- 



ization model to the thermal component (Townsley et al. 
2006). The thermal component is found to be neon- rich 
and is also slightly over-abundant in oxygen and magne- 
sium, in respect to LMC values, indicating that much of 
the thermal emission arises from the SN ejecta. 



To see where the metal-rich materials are located, we 
present the EW images of the O, Ne, and Mg lines in 
Fig. These images are constructed in a method sim- 
ilar to those used by Hwang et al. (2000) and Park et 
al. (2002). A major difference of our method from theirs 
is that we rebin the data using an adaptive mesh with 
each bin including at least 10 counts(see Warren et al. 
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FlG. 4. — Raw ACIS-S image of N157B. The labeled regions 
are used for spectrum extraction. The small circle in Region PI 
encircles the emission of PSR J0537— 6910, and other point sources 
are removed. The coordinate axes are for R.A. (J2000) and Dec. 
(J2000), respectively. 




energy (keV) 

Fig. 5.— The overall spectrum of SNR N157B fitted with the 
model of a power-law plus an NEI thermal plasma. The lower panel 
(likewise also in other spectral plots) shows the residuals of the data 
minus the model in units of cr. The two spectral components are 
also shown separately. 



2003 for similar binning). The figure shows that the Une 
emission is distributed over a broad region, except for the 
cometary X-ray nebula, which is dominated by nonther- 
mal emission and appears featureless. The Ne emission 
(in green) has the largest extent, while the Mg line (in 
blue) is strong mostly near the southwestern edge of the 
nebula (Region Dl), while the O line (in red) is concen- 



TABLE 3 
Spectral Fit for the Entire SNR 



parameter 



value 



xVd-o.f. 

Na (lO^l cm-2) 



299.2/301 
6.1 ± 0.6 



power-law component 



flux (10 ^^crgscm ^s ^) 

NEI thermal component 

kT (keV) 

[O/H] 

[Ne/H] 

[Mg/H] 

[Si/H] 

n^ti (10l°cm~3s) 
fnenaV^ (10^** cm'^) 
flux (10""'^"'^ ergscm"^ 



2.29 
1.4 



+0.05 
0.06 



0.72-' 



,+0.13 
-0.10 

43+0 16 

1 91+0.44 
^•^^-0.13 
0^7g+0.37 



-0.19 
-,+0.40 
-0 
a+3.05 
-1.16 
-1.49 
-1.95 
0.7 



0: 

2.39:! 
5.22"* 



Note. — The unabsorbed fluxes are in the 0.5-10 keV band. 
The net count rate is 1.090 ± 0.005 counts s~i. 

^ Quantity fneUnV is the volume emission measure of the hot 
gas, where / denotes the filling factor. 



trated in a bright northeast patch (D3). 



0.01 



E \ \ \ \ 











0.5 



energy (keV) 



Fig. 6.— ACIS-S spectrum of PSR J0537-6910 and the best-fit 
power-law. 



2.2. PSR J0537-6910 

The location of PSR J0537— 6910 corresponds to the 
emission peak of the N157B nebula. The position listed 
in Table El is consistent with the time-resolved peak 
of the pulsed emission from the pulsar (5''37™47^.36, 
69°10'20".4; Wang et al. 2001; see also Mignani et al. 
2005). We extract a spectrum of PSR J0537-6910 or 
the core of the nebula bar from a circle of 1" radius (90% 
EER) centered on the emission peak and a background 
from the bar minus the circle (Region PI in Fig. ^J. 
We fit the spectrum with a power-law model (Fig. 
no evidence is found for any significant pile-up, based 
on the model given by Davis 2001). The fit is satisfac- 
tory (x^/d.o.f. = 110.0/125) and gives the best-fit power 
law index as F = 1.73toJg ^^"^ foreground absorp- 



tion column density as 



5.6ij;:^ X 10" 



The 
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absorption-corrected flux is ^ 1.9 x 10^^^ ergscm^^ s^^ 
in the 0.5-10 keV band (or a luminosity of ~ 5.7 x 
10'^^ ergs s""'^). The results are consistent with the pre- 
liminary results given by Mignani et al. (2005). The 
power law index is also consistent with that measured for 
the pulsed X-ray emission from the pulsar (~ I.G^q'j), 
based on the ASCA and RXTE spectra (Marshall et al. 
1998). Assuming the same power law (F = 1.73), we find 
that the pulsed fraction is 47% in the 0.5-10 keV band, 
compared to 40% in the 0.2-4 keV band (Wang et al. 
2001). 

In the Spitzer IRAC images there are no point sources 
at the location of the pulsar. The upper limits for the 
detection of a point source with 3a confidence are flux 
densities of 22, 28, 53, and 143 ^Jy in the IRAC 3.6, 4.5, 
5.8, and 8.0 /im bands, respectively. Similarly, no radio 
or optical counterpart has been detected for the pulsar 
(Crawford et al. 2005; Mignani et al. 2005). 

2.3. Diffuse Emission: Nonthermal Component 

Our analysis of the global diffuse X-ray distribution 
shows that the cometary nebula has a hard and nearly 
featureless spectral characteristics, confirming the con- 
jecture made in the previous studies (e.g., WG98). Fig.|3 
further shows some interesting details that have not 
clearly been identified before. These include evidence 
of a halo around the bar that is most apparent in front 
of the bar and extends about 5" from the pulsar. In 
the 0.3-1.5 keV band, the new data also show consid- 
erable substructures in the tail. However, the counting 
statistics of the data does not allow for a quantitative 
characterization of these substructures. 

We examine the spectral properties of the key non- 
thermal features, individually. Fig. |H1 shows the net 
spectrum of the halo, extracted from Region PO (a lo- 
cal background described in § 12.41 is subtracted). In 
addition to a dominant contribution from a power law 
with a photon index of r~ 1.7 — 2.1, there is clearly a 
thermal component, as evidenced by the presence of the 
~ 1.34 keV Mg XI Ka line. Characterized by a thermal 
plasma {kT 0.6-1.2 keV), this component accounts for 
about 20% of the total emission from Region PO. The 
exact contribution from this thermal component is un- 
certain, somewhat depending on the background sub- 
traction. The component, rich in magnesium (~ 0.6-4.5 
solar), may represent an enhanced surrounding thermal 
plasma emission projected along the line of sight. 

Fig. ini shows the net spectra of regions PI (with the 
pulsar excluded) and P2-P8 (Fig. gj). The PI spectrum 
is dominated by the SW-NE oriented bar (- 2.4" x 6.8"), 
which has a very high-surface brightness. We can thus 
divide the bar region into individual segments, but find 
no statistically significant spectral variation on either 
side of the pulsar. The spectrum of Region PI can be 
well fitted (x^/d.o.f. = 206.3/191) by a power-law with 
a photon index of 2.18^0 53 and a foreground absorp- 
tion Nil ^ 6.1^2^4 X lO'^^cm"^. The accumulated spec- 
trum of the P2-P8 regions may also be characterized 
approximately (x^/d.o.f. = 318.1/241) by a power-law 
with an index of 2.62 ± 0.04 and a hydrogen column 
JVh ~ (6.3 ± 0.2) X 10^1 cm-'^. The N^ values are consis- 
tent with those inferred from the spectral analyses of the 
entire remnant I2.1|l and the pulsar fS I2.2|l . The spec- 



tra from the individual P1-P8 regions, fitted well with 
the power law, show a systematic change of the photon 
index along the major axis of the nebula fFig. llO|l . Here 
we have fixed N^ to the value 6.1 x 10^^ cm^^; if N^ is 
allowed to vary in the fit, the result is very similar, albeit 
with larger error bars. 



2.4. Diffuse Emission: Thermal Component 

We characterize the properties of the diffuse X-ray 
emission away from the nonthermal cometary nebula by 
analyzing the ACIS-S spectra of individual regions of 
N157B (Fig. m Fig.ini)- The spectrum of Regions Ll- 
3, excluding the background estimated in Regions Bl 
and B2, is used to characterize the general low-surface 
brightness emission from the SNR. For the spectral anal- 
ysis of the substructures (Dl, D2, and D3), we adopt 
the Ll-3 spectrum as the local background. All these 
spectra show distinct emission lines, indicating predom- 
inant contributions from optically-thin thermal plasma. 
Indeed, the spectra can all be fitted well with the NEI 
model. Similar to the analysis of the accumulated spec- 
trum from the entire SNR f^i I2.1|l . the model fits allow 
the abundances of key line elements to vary, while leaving 
the rest of the metal abundances fixed to the LMC value 
(0.3 solar). Table 01 summarizes the net count rates and 
results of the spectral fits, as well as the inferred energy 
fiuxes and hydrogen number densities of the individual 
regions. The densities are derived on the assumption 
that the 3-D shapes of the elliptical regions (D1-D3) are 
oblate ellipsoids and those of the rectangular regions (Ll- 
L3) are cuboids with the average line-of-sight size to be 
the remnant radius. The deviation from this assumption 
and the non-uniformity of the X-ray-emitting plasma are 
consolidated into the "filling factor," /, of the individual 
regions. The regions show significant spectral variation 
among them. While the metal abundances in Regions 
Ll-3 are consistent with the LMC value, the substruc- 
tures all exhibit significant metal enrichment and vari- 
ations in physical conditions. For example. Region D3 
shows both a strong O VIII Ka line and a low plasma 
temperature, which explains the softness of its emission 
(e.g.. Fig. n^)- Regions Dl and D2 show comparable 
temperatures, densities, and neon abundances, but Re- 
gion Dl is also particularly abundant in Mg, seen in both 
the spectrum and the EW map, and its ionization time- 
scale is significantly smaller than that of Region D2. 

It is possible that part of the diffuse X-ray emission 
may have origins other than the N157B SNR. However, 
there is little evidence for a spatial correlation of the 
diffuse X-ray intensity with stellar concentrations (e.g.. 
Fig. ^), indicating that the direct stellar emission and 
the gas heating due to stellar winds from massive stars 
(e.g., wind- wind and wind-cloud collisions) do not con- 
tribute significantly. The X-ray emission from the pre- 
existing superbubble is substantially softer than what is 
inferred from the SNR and has a very low-surface bright- 
ness (Smith & Wang 2004) . Our analysis is based on the 
enhanced emission above the local background, which in- 
cludes the superbubble contribution. 

3. DISCUSSION 

Based on our analysis, we now attempt to present a 
coherent interpretation of the characteristics of N157B: 
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RIGHT ASCENSION (J2000) RIGHT ASCENSION (J2000) 

Fig. 7. — Close-ups of the PWN in the 0.3-1.5 keV band (a) and the 1.5-7 keV band (b). These raw count distributions, smoothed with 
a Gaussian with FWHM=0''3, are plotted logarithmically. 



TABLE 4 

Spectral Fits for Individual Thermal Features 



parameter 


Ll-3 


Dl 


D2 


D3 


net count rate (10~^ counts s~^) 


1.79 ± 0.10 


4.06 ±0.11 


3.70 ±0.10 


3.76 ±0.13 


xVd.o.f. 


23.6/24 


84.3/60 


87.1/62 


45.4/40 


Nh (10^1 cm-2) 


3.8 ±0.6 


7 4+1-0 
'•^-0.4 


4.81°:? 


9-3+28 


kT (keV) 


1 04+°-*^ 
^■"^-0.31 


1 47+0.59 


1 99 + 0.41 
J^-^^-0.31 


23+0-15 
'-'•^'-'-0.07 


[O/H] 










[Ne/H] 




1 OO+0.36 
r.JO.Q.so 


0.96 ±0.26 


'Sl+O-50 


[Mg/H] 




r, r,q + 1.10 

^-^'-'-0.84 


0.501°:?^ 




neti (10^°cm-3s) 


c oq+8.30 


0-68t°;f5 


9 oQ+0.93 
^••JM_0.96 


2.30i»:l| 


frieriiiV (10^7 cm"^) 


3.20ti'j'5 


7 oq+4.42 
' -■^'-'-2.50 


"'•'0-1.44 


540{< 3460) 


flux (10~^^ ergscm"-^ s~^) 


2.0 


14.5 


6.0 


313 


na/l/^ (cm-3)'' 


0.22 


0.56 


0.67 


3.8 



Note. — The fluxes are in the 0.5-10 keV band. 
^ In the estimate of the densities, we assume cuboidal volumes for rectangular regions LI— L3, with sizes 8.0 X 7.5 X 12 pc'^, 14 X 5.2 X 12 pc^, 
and 7.3 X 2.9 X 12 pc^, respectively, and oblate spheroids for elliptical regions D1-D3, with half axes 8.6 X 8.6 X 2.2 pc^, 5.9 X 5.9 X 2.1 pc^, 
and 7.2 X 7.2 X 5.0pc3, respectively. 



the exceptionally large sizes of both the thermal and non- 
thermal components and the lack of a well-defined outer 
shell as well as the distinct spatial/spectral substructure 
and metal enrichment pattern. 

3.1. Environmental effect 

Nf 57B is presumably associated with the OB associa- 
tion LH99, as they have comparable foreground absorp- 
tion/extinction (WG98), in addition to their apparent 
superposition in the sky. While the characteristics of 
N157B appear to be a manifestation of its association 
with LH99, the new X-ray results provide new insights 
into the nature of the SN progenitor and the proper- 
ties of the ambient environment. First, we can use the 
observed pattern of metal enrichment to constrain the 
mass of the progenitor star. As the presence of a pul- 



sar implies a core-collapse nature of the SN, we plot 
in Fig. 1121 the relative metal abundances predicted from 
such SN models for stars with low metallicities (Thiele- 
mann, Nomoto, & Hashimoto 1996, hereafter TNH96; 
Woosley & Weaver 1995, hereafter WW95). The enrich- 
ment is compatible with the TNH96 model for a pro- 
genitor mass of 2O-25M0, or the WW95 models for a 
progenitor mass of 30-35Mq, Z — O.IZq, and an ex- 
plosion energy of ~ 1.2 x 10^^ ergs. Although these two 
model predictions do not overlap, they both indicate a 
progenitor mass of ^ 20Mq. Recent models taking into 
account stellar mass loss suggest that neutron stars can 
be produced by metal-poor progenitors only if the pro- 
genitor mass is < 25Mq (Heger et al. 2003). This mass 
limit, while consistent with the above estimate from the 
TNH96 model, is lower than the masses of the 03 and 05 
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Fig. 8. — ACIS-S spectrum in the regions PO, fitted with a 
power-law plus an NEI thermal plasma. 
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Fig. 9.— ACIS-S spectra in the Regions PI (a) and P2-P8 (b), 
together with power-law fits. 



Stars that still exist in LH99 (Schild & Testor 1992). The 
earlier explosion of a less massive star thus implies that 
the star formation in the vicinity of LH99 is not coeval. 
Further evidence for a prolonged star formation has re- 
cently been provided by the massive proto-stars revealed 
in Spitzer observations; for example, the brightest point 
source in the Spitzer 8 /im image is a proto-star (Fig.^). 
We next examine how the co-existence with the OB 




10 20 30 40 50 



z (arcsec) 

Fig. 10. — Variation of the power-law photon index of the non- 
thermal tail as a function of the distance from the pulsar. Curves 
are plotted according to the particle ejection model represented by 
Eq. {TOJ with various u and Bj^ parameter combinations: 0.5c and 
0.9 X 10— ' G (dotted curve), 0.6c and 0.9 X 10"'' G (short-dashed), 
0.6c and 1.0 X 10"* G (sohd), 0.6c and 1.1 X 10"" G (long-dashed), 
and 0.7c and 1.1 X 10""' G (dot-dashed). 



association LH99 may have influenced the evolution of 
N157B. Fig.[T]i shows a Spitzer IRAC 8 /im view of the 
large-scale environment of N157B. The 8 /im band image 
is dominated by emission from polycyclic aromatic hy- 
drocarbons (PAHs). As PAHs can be easily dissociated 
by UV radiation from O stars, the 8 /im emission traces 
H I gas just outside the ionization front. The 8 /im im- 
age shows that N157B is inside a large cavity that might 
be a superbubble produced by the mechanical energy in- 
put from the OB association LH99. The existence of 
a blister-type superbubble around N157B has also been 
suggested previously based on the velocity structure of 
dense ionized gas (Chu et al. 1992). This IR cavity might 
also be connected with the cavities and "tunnels" pow- 
ered by multiple OB associations, including also LH90 
to the east, as indicated by the presence of a large-scale 
(~ 20') low-surface brightness X-ray emission in the re- 
gion (Smith & Wang 2004). The exact shape of the mid- 
infrared cavity, let alone its boundaries, is very uncertain; 
there may be substantial projection effect. In morphol- 
ogy, the pulsar, presumably close to the explosion center, 
is located at the west edge of the cavity. The large-scale 
diffuse X-ray emission extends into the cavity. The ex- 
pansion of the SNR into such a low-density hot medium 
at least partly explains its unusual large size and its lack 
of a well-defined outer shell. The hot ambient medium 
has a high isothermal sound velocity, thus the shocks are 
expected to be weak with small Mach numbers, and the 
blastwave may become thermalized and deposit the bulk 
of the SNR energy into the thermal energy of the hot 
medium (Tang & Wang 2005). The relatively low pres- 
sure inside the SNR further allows for the expansion of 
the PWN to an exceptionally large volume. 

While the overall morphology of N157B is apparently 
affected by the large-scale superbubble environment, the 
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energy (keV) energy (keV) 

Fig. 11. — ACIS-S spectra of selected thermal substructures (see Fig. Bl LI— 3, Dl, D2, and D3). All these spectra are fitted with the 
NEI model (Table and are re-binned to achieve a background-subtracted S/N ratio of 3. 



presence of a dense cloud close to the explosion site seems 
to play a critical role in shaping the dynamics of the ma- 
terials around the pulsar. This provides an alternative 
explanation for the trailing PWN structure, which we 
previously proposed to be due to pulsar motion (WG98; 
Wang et al. 2001). Probably, the trailing tail is not 
necessarily produced by the passage of the moving pul- 
sar; and an existing example may be the Duck X-ray 
Nebula, where the tail is ejection of pulsar wind par- 
ticles at weak relativistic bulk velocity confined by the 
ambient medium ram-pressure (Kaspi et al. 2001). In 
SNR N157B, the pulsar appears to be near the explo- 
sion center and have not moved out of the SNR, at least 
in projection. The size of the nonthermal nebula, es- 
pecially its head, is large, which is somewhat difficult 
to explain in a supersonic relative motion scenario, al- 
though the physics of the shocked pulsar wind is quite 
uncertain, complicated by its non-spherical symmetry. 
Instead, the subsonic ram-pressure confinement can be 
consistent with the gas motion relative to the pulsar (see 
below and WS.2\i . The cloud is partly seen as the strongest 
diffuse mid-infrared emission peak south to the pulsar in 
Fig-Gt- The cometary PWN tail points to the northwest, 
almost parallel to the long dimension of the cloud or the 
southwest edge of the superwind cavity (Fig. la), con- 
sistent with the relatively narrow extent of the thermal 
emission on the southwest side of the PWN (Fig. lb and 



Fig. Id). The northern part of the cloud, not apparent in 
Fig.^, is bright in optical (Fig.^). The intensity of the 
enhanced diffuse optical light shows a reasonably good 
morphological similarity to the relatively high surface- 
brightness X-ray emission in the 0.3-1.5 keV band. The 
dust lane at the southeast side of Fig. ^ casts a shadow 
in the emission distribution and corresponds to harder 
X-ray emission (Fig. Qd), apparently due to the X-ray 
absorption (but the limited counting statistics and the 
uncertainty in the background subtraction make it hard 
to quantify the absorbing gas column). In fact, the as- 
sociation of N157B with this optical nebula has already 
been suggested by Chu et al. (1992), based on its kine- 
matics. In view of these observational factors, we suggest 
a physical association that is produced by the impact of 
N157B on the cloud. Although the shape and boundary 
of the mid-infrared cavity are very uncertain (consider- 
ing substantial projection effect) and the exact interplay 
(the structure and dynamics) between the two cannot 
yet be determined, this impact likely led to a strong re- 
flection shock, which has probably swept through and 
been further reflected around much of the PWN. It is 
also possible that there is a reflection from the east edge 
of the cavity, although we find little optical evidence 
for a strong interaction. The resultant heating by re- 
flection shock and possibly large bulk flow of the SNR 
gas may thus be responsible for both the observed one- 
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sided morphology of the PWN and the thermal emission 
enhancement (like a sheath) around the southern part 
of the PWN. The enhancement east of the pulsar (Re- 
gion D2) is particularly strong just (^ 10" — 20") east of 
the pulsar and is bordered by a dust lane farther to the 
east. It shows no evidence for nonthermal emission. This 
part of the enhancement is clearly edge-enhanced and 
has corresponding optical brightening (Fig. ^), proba- 
bly representing the strongest ongoing reflection of the 
SNR from the cloud. The high-pressure of the reflected 
gas may have greatly influenced the confinement of the 
pulsar wind and hence the direction of its outflow toward 
the northwest. This interpretation of the PWN and its 
surrounding medium is similar to the scenarios proposed 
by WG98 and van der Swaluw et al. (2003): the observed 
one-sided morphology is a result of an asymmetric ram- 
pressure confinement. But, the present interpretation 
of the strongly one-sided PWN morphology does not re- 
quire a fast proper motion of the pulsar; the ram-pressure 
(or the relative motion between the pulsar and the sur- 
rounding medium) can naturally arises from the highly 
asymmetric reflection shock of the SNR (see more dis- 
cussion below). Having proposed this alternative, we do 
not abandon the pulsar proper motion scenario, which is 
still plausible. 

3.2. Evolutionary Stage of SNR N157B 

The properties of N157B is determined not only by 
its environment but also by its specific evolutionary 
stage, although its age is still very uncertain. The spin- 
down age of the associated pulsar PSR J0537— 6910 is 
~ 5 X lO'^yr, as derived from the pulse period and 
its derivative, assuming a negligible initial pulsar pe- 
riod and a dipole magnetic field configuration (Mar- 
shall et al. 1998). These assumptions may not be valid 
for young pulsars in general and for this pulsar in par- 
ticular because it displays repeated episodes of strong 
timing glitches (Marshall et al. 2004). Nevertheless, 
the spin-down age is roughly consistent with the large 
size of the remnant (12 pc in radius for a total angu- 
lar extent of ^ 100"). From the ionization timescale 
{rieti ~ 2.4 X 10^° cm^'^ s; TableOl and the mean gas den- 
sity (nn - 0.45/~i/2j^j"2^/2gjn-3 ^Yieie R12 = R/12pc, 
assuming rig ^ 1.2nH), we may estimate an average ion- 
ization age of the X-ray-emitting plasma in the rem- 
nant (i.e., the time elapsed since it was shocked) as 

ti ^ 2 X 10'^/^/^i?]^2^ yr, which is considerably shorter 
than the pulsar spin-down age, especially because / can 
be small. This discrepancy may be understood if the X- 
ray emission originates primarily from the plasma that 
is recently shock-heated. 

Indeed, the thermal properties of N157B, as inferred 
from § 2.4, suggest that the X-ray-emitting gas consists 
of various metal-enriched clumps embedded in a more 
diffuse plasma. This diffuse plasma is distributed more 
widely than the clumps and also seems to have both a 
nominal LMC metal abundance and a larger average ion- 
ization age, which is comparable to the spin-down age of 
the pulsar. Therefore, the plasma may mostly represent 
the shocked ISM, whereas the clumps are mainly due to 
SN ejecta. Unfortunately, the quality of the data is not 
adequate to make a quantitative decomposition of these 
two components. We could only make a crude estimate 



of the thermal energy as ~ 3 x 10^°/"'^/^i?ij2^ ergs and the 

mass as ^ 110/^/^i?i^2^MQ. This estimate is subject to 
the uncertainty in the effective filling factor, which can- 
not be determined and is small, probably on the order 
of ^0.1. This gas mass is comparable to several solar 
masses expected for the SN ejecta and is thus consistent 
with the significant metal-enrichment detected. Depend- 
ing on the / value, the thermal energy may be consider- 
ably smaller than the canonical SN energy (~ 10^^ ergs). 
A substantial fraction of the energy may still be in the 
kinetic form, carried primarily by a weak outward blast- 
wave into the low-density hot superbubble interior (§ 3.1; 
Tang & Wang 2005). 

Among the substructures considered, the enhance- 
ment in Region D3 shows distinct characteristics: a par- 
tially round shape with a size of ~ 30" or 7pc, a dis- 
tinctly low gas temperature (~ 0.2 keV), a high density 
(~ 4/~^/^cm~^), a low ionization age (~ 160/^/^yr), 
and an apparent overabundance in oxygen. This fea- 
ture may represent a clump of oxygen-rich ejecta recently 
heated by a low- velocity reverse shock 400 km s"'^), 
as inferred from the gas temperature. However, this hy- 
pothesis cannot explain the morphology or size, as the 
ionization time scale (Table 4) would suggest a velocity 
of ~ 4 X 10**/^^/^ km s^^, much higher than the value 
inferred from the plasma temperature. Of course, the 
NEI plasma we used may be an over-simplification for 
the modeling of the spectrum fFig.lll|. and the inferred 
ionization time scale may be misleading. Without the 
constraint on the time scale, the D3 enhancement may 
simply be part of the shocked SN ejecta, possibly bal- 
looning out from the high-pressure SNR interior into the 
superbubble. In this case, the relatively low tempera- 
ture can be naturally explained by the adiabatic cooling 
of the expanding gas. Alternatively, D3 may arise from a 
separate oxygen-rich SN, which is quite possible with the 
presence of the OB association LH99. In this case, the 
partially round morphology can naturally be explained as 
a coherent shell of reverse-shocked SN ejecta. The age of 
the remnant should be comparable to the ionization age 
inferred above. The reverse shock velocity is expected 
to be substantially lower than that of the forward blast- 
wave, which is probably very weak in X-ray emission. If 
the remnant is also expanding inside the superbubble, 
the lack of an apparent radio or optical counterpart is 
then expected. The absence of a PWN associated with 
the remnant further suggests that the stellar progenitor 
of the SN is more massive (> 25Mq; Heger et al. 2003) 
and that the stellar remnant may be a black hole. In 
fact, we are also a bit uncomfortable about involving a 
second SNR in the interpretation. But if we accept the 
spectral results (mainly the small ionization parameter) 
from the NEI model fit to the D3 spectrum, we find no 
better solution than a new SNR. 

As discussed in § 3.1, a "sheath" of the thermal X- 
ray enhancement (including D1-I-D2) around the PWN is 
likely a result of its interaction with the material reflected 
off the dense cloud. The physical processes involved are 
similar to those in the interaction between an SNR re- 
verse shock and a PWN of a fast moving pulsar, as sim- 
ulated by van der Swaluw et al. (2004). Before such an 
interaction, the PWN consists of two main components, 
the shocked pulsar wind material and the swept-up SN 
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Fig. 12. — The best-fit abundances relative to Mg, compared 
with the predictions from core-collapsed SN models (TNH96 and 
WW95). Note that because the best-fit Si value is only ~ 10~*, it 
is invisible in the right panel (with a logarithmic scale). 



ejecta between the forward shock and the contact discon- 
tmuity (Reynolds & Chevaher 1984). The thermal emis- 
sion from the swept-up SN ejecta is expected to be weak 
and is confined within a relatively thin shell. The pas- 
sage of the strong reflection shock (e.g., stage 2 in van der 
Swaluw et al. 2004) compresses the PWN, sweeps part of 
it into the trailing tail, and is meanwhile further reflected 
forward. Consequently, both the previous undisturbed 
and the swept-up SN ejecta are compressed and heated at 
the front side of the PWN. The density and temperature 
in this "sheath" around the PWN should be enhanced 
by a factor of 1-2.5 (compared to those in the upstream 
gas), depending on the Mach number of the reflection 
shock relative to the pulsar (Spitzer 1982; Hester & Cox 
1986). Such density and temperature contrasts are con- 
sistent with our measurements, albeit with large uncer- 
tainties (TableEJ. We thus conclude that the interaction 
between the reflected material and the PWN provides 
a natural interpretation of both its strongly one-sided 
morphology and its surrounding thermal X-ray emission 
enhancement. 

3.3. Physical Properties of the PWN 

Based on the spectroscopic results of various nonther- 
mal components of N157B, we attempt to determine their 
nature and quantify their physical parameters. The bar 
around the pulsar PSR J0537— 6910 has been proposed 
to be a torus (Wang et al. 2001), similar to those de- 
tected around several well-known young pulsars [e.g.. 
Crab (Weisskopf et al. 2000), Vela (Helfand, Gotthelf, 
& Halpern 2001), PSR 0540-69 (Gotthelf & Wang 2000), 
and G54. 1-1-0.3 (Lu et al. 2002)] and beheved to repre- 
sent the termination shocks of the toroidal pulsar winds 
(e.g., Ng & Romani 2004). This proposal is consistent 
with the results from our spectral analysis of the bar. 
The lack of a systematic spectral variation across the 
bar simply reflects similar particle acceleration and syn- 
chrotron emission throughout the termination shock re- 
gion, just as in the Crab nebula (Mori et al. 2004). The 



power-law photon index inferred from the analysis also 
agrees well with the values (~ 2) observed in the other 
toruses (Gotthelf 2003). The centrally enhanced bar- like 
morphology can be naturally explained by the edge-on 
orientation of the torus because of the expected Doppler 
boosting, which depends on the angle between the bulk 
motion velocity of shocked pulsar particles and our sight 
line (e.g., Komissarov & Lyubarsky 2003). Similar ef- 
fects have been observed in both the Crab nebula and 
G54.H-0.3 (Lu et al. 2002). The size of the torus should 
be comparable to the radius of the termination shock 
(modifled from Chevalier 2000), 




(1) 



where e b denotes the magnetic field fraction of the pulsar 
wind luminosity (~ 0.5 in the equipartition case) and Ot 
is the half angular span of the toroidal wind (i.e., 47r sin 9t 
is the solid angle subtended at the pulsar by the torus). 
The pulsar spin-down power E = 4.8 x 10'^* ergs s~^ 
(Wang et al. 2001; Mignani et al. 2005) is adopted. The 
magnetic field B just behind the shock is not well un- 
derstood. But if it is comparable to the average strength 
(- 2 X 10"-* G) in the PWN and 0t = 20° (as estimated 
in Wang et al. 2001), we can then estimate r* w 0.6 pc or 
2'.'4. This estimate is comparable to the observed overall 
dimension (^ 2'.'4 x 6'.'8) of the bar, and thus both the 
morphological and spectral characteristics of the bar are 
consistent with its being an edge-on torus resulting from 
the terminal shock of the pulsar toroidal wind. 

The N157B PWN as a whole shows a featureless and 
hard spectrum. The spectral steepening along the ma- 
jor axis of the tail indicates significant synchrotron cool- 
ing. Similar spectral properties have been seen in other 
SNRs, such as G21. 5-0.9 (Slane et al.. 2000; Warwick 
et al. 2001, Safi-Harb et al. 2001), 3C58 (Torii et al. 
2000), IC443 (Bocchino & Bykov 2001), and possibly 
in W44 (Petre, Kuntz, & Shelton 2002); however, their 
X-ray softening has not been quantitatively compared 
with models. The bright and long trailing pattern of the 
N157B nebula makes it possible to quantitatively com- 
pare observations with models. 

We construct a simple 1-D toy model, similar to 
the spherically symmetrical treatment by Amato et al. 
(2000). In the 1-D approximation of the N157B PWN, 
the shocked wind particles stream off the bar along the 
downstream tail at a speed u within a cylinder of radius 
r. The differential particle density distribution as a func- 
tion of the distance from the injection point (the bar) z 
is then 



NiE,z) 



J(Sini) 



dt. 



dz 



dE 



(2) 



where E'ini and tini are the initial energy and time of 
the injected particle, while the injection rate J(£'ini) is 
assumed to be a power law, 

J(i?i„i) = KErP. (3) 

If the bulk velocity u is assumed to be constant (i.e., 
z/t), the energy loss of the particles resulting from the 
synchrotron radiation is then 



dE 
'dt 



-ciBlE\ 



(4) 
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where ci = 4e^/(9rngC^) and B± is the magnetic field 
component perpendicular to the particle motion. Inte- 
gration of the above equation gives E^ini — E / {\ — E / Eoi), 
where 

E^ = u/ic^Blz) (5) 

marking the maximum energy of the particles after the 
synchrotron loss at z. Then the particle density distri- 
bution is given by 



N{E,z) 



K 



E-P 1 



E 

Em 



p-2 



(6) 



The synchrotron emissivity of the particles is je = 
{ciB\E'^)N [E , z)dE/de, where the photon energy is re- 
lated to the particle energy as e = C2hB^E^ [where 
C2 = 0.29 X 3e/(47rm3c5)]. Therefore 



Cl 



K 



2c2h: 



BiE^~P 1 



E 

Em 



p-2 



(7) 



The dependence of the emission's photon index, 

r = -(1 + dlnje/dlne), (8) 

implicitly on e and z (via E and Em, respectively) is 
derived as 



r 



p+l E/En 



2 1 - E/En 



(9) 



The above equation shows that the spectral index varies 
with both the photon energy and the distance from injec- 
tion point. Since the spectral indices are determined in a 
broad energy range, to compare with the observed values, 
we calculate a weighted spectral index, f = jmmCeu^^m) 
[jij^, ^)/e]r(i', z) de/ J[j{iy, z)/e] de, namely 

p+l p-1 (V^m/ei - l) - (^v^em/niin(£u,em) 



ei-1 



p-i 



(10) 

where ei and £u are the lower and upper energy limits of 
the spectrum used, while £m is the photon energy corre- 
sponding to Em, 



17 pc 



B, 



10-4 G 



keV. 



(11) 



Eq. 1101 can be compared with the measured values. 
We fix the photon index at the injection point (PI), 
r = 2.17 (for iVn = 6.1 X 10^^ cm~^), which corresponds 
to p = 3.34, and adopt ei = 0.5 keV and Eu — 7keV. 
Fig. 1101 shows predicted curves, assuming the magnetic 
field to be around the mean PWN value (Wang et al. 
2001). These curves match the measured values well, 
suggesting u^0.5c, comparable to the sound speed of 
an ultra-relativistically hot plasma (cs = c/^/3). Simi- 
lar high downstream flow speeds are observed or inferred 
at the terminal shocks of the Crab nebula (Hester et 
al. 2002) and G54.1+0.3 (Lu et al. 2002); but in these 
cases, the flows are presumably quickly damped and en- 
ergy is then transported primarily due to particle diffu- 
sion. In an one-sided pressure confined PWN, the down- 
stream fiow is re-directed into a bulk motion at a sim- 
ilar or higher speed, as demonstrated analytically (e.g., 



Wang, Li, & Begelman 1993) and in simulations (e.g., 
Bucciantini et al. 2005). Therefore, a strong pressure- 
confinement of the PWN is needed to give a consistent 
interpretation of both the exceptionally large linear size 
and the observed spectral steepening of the nonthermal 
component of the N157B PWN. 

4. SUMMARY 

Based on our on-axis Chandra ACIS observation, we 
have conducted a detailed spatially-resolved X-ray spec- 
troscopy of N157B. Combining the X-ray analysis with 
the optical and mid- infrared observations, we have ex- 
plored the physical origins of various distinct character- 
istics and their implications. Our main results and con- 
clusions are summarized as follows: 

1. We detect eight point-like X-ray sources in the 
ACIS observation. Two of these sources in the 
vicinity of the pulsar PSR J0537— 6910 have near- 
IR and optical counterparts, which are identified 
as massive stars in the LMC. The relatively high 
X-ray luminosities of the two sources indicate that 
they are probably colliding stellar wind binaries. 
A comparison of the X-ray and 2MASS positions 
of the two sources suggests that the absolute as- 
trometry of the ACIS observation, hence the X-ray 
centroid position of the pulsar, is better than about 
0'.'5. 



2. We confirm the nonthermal nature of the various 
high-surface brightness X-ray features, such as the 
bar, tail, and halo, as suggested in the previous 
morphological and integrated spectral studies. The 
bar shows no significant spectral shape variation 
and most likely represents the terminal shock of a 
toroidal pulsar wind. The strong spectral steepen- 

p-2 ing observed in the tail is consistent with the syn- 
chrotron radiation loss and suggests a bulk velocity 
p-i of outfiowing pulsar wind particles to be ^ 0.5c. 

3. The thermal emission is unambiguously resolved 
to arise from a region much larger than the non- 
thermal PWN and accounts for about one-third of 
the total X-ray emission in the 0.5-10 keV band. 
The thermal emission consists of various clumps 
embedded in a diffuse X-ray-emitting plasma. The 
clumps are over- abundant in a elements and thus 
represent the SN ejecta. The over-abundance pat- 
tern suggests a mass <; 20AfQ for the SN progeni- 
tor. However, a partially-round-shaped soft X-ray 
clump in the outer region of N157B shows dis- 
tinct spectral characteristics, unusually low in both 
temperature and ionization time scale and over- 
abundant in oxygen. This clump may thus rep- 
resent a separate young SNR. 

4. The thermal emission of N157B is enhanced around 
the nonthermal PWN, particularly near its head. 
This enhancement is perhaps a result of the recent 
compression and heating of the ejecta due to the 
SNR reffection from a nearby dense cloud. This in- 
teraction of the SNR reffection and the PWN nat- 
urally explains its one-sided morphology. 

5. The relatively weak thermal emission, the lack of a 
rim-brightened outer blastwave, and the large sizes 
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of both the thermal and nonthermal components all 
point to a scenario that SNR N157B is expanding 
into a low-density hot medium, which likely repre- 
sents the interior of a superbubble created by the 
co-existing OB association. 
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